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Abstract 

We have found the graviton contribution to the one-loop quantum correction to the 
Newton law. This correction results in interaction decreasing with distance as 1/r 3 and 
is dominated numerically by the graviton contribution. The previous calculations of this 
contribution to the discussed effect are demonstrated to be incorrect. 
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1 Introduction 



The problem of corrections to the equations of motion, arising in general relativity, is far 
from being new. The classical relativistic corrections to these equations were found long ago 
by Einstein, Infeld and Hoffmann an d by Eddington and Clark 0. (A relatively simple 
derivation of these corrections is presented in the textbook jj.) Later this result was reproduced 
by Iwasaki by means of Feynman diagrams @]. Thus, the problem of the classical relativistic 
corrections to the Newton law is solved finally []. 

Let us note that the general structure of the relativistic classical correction to the interaction 
potential of two bodies with masses mi and m.2, which would be of second order in the Newton 
gravitational constant k, is clear immediately. Indeed, the quantity km/c 2 (c is the velocity 
of light) has the dimension of length, so that with the account for the symmetry under the 
interchange mi <-> vri2 the correction should be of the form 

TT k 2 m 1 m 2 (m 1 + m 2 ) , lN 
U <* = a d ^2 • W 

The dimensionless constant a c \ as found in the above works equals 1/2. 

There is one more linear in k combination of constants which can be used for the construction 
of a power correction to the Newton potential. We mean 

kh 2 

where h is the Planck constant, l p = 1.6 • 10~ 33 cm is the Planck length. Clearly, such a 
correction, being of course of quantum nature, should look as follows: 

_ k 2 hm l m 2 
u qu - a qu — . {l) 

One has to find the numerical constant a qu . In spite of extreme smallness of the quantum 
correction, its investigation certainly has a methodological interest: this is a closed calculation 
of a higher order effect in the nonrenormalizable quantum gravity. 

The reason why this problem allows for a closed solution is as follows. The Fourier-transform 
of 1/r 3 is 

V Xp( 7 qr) =-2,\W- (3) 

/yd 

This singularity in the momentum transfer q means that the correction discussed can be gen- 
erated only by diagrams with two massless particles in the t-channel. The number of such 
diagrams of second order in k is finite, and their logarithmic part in q 2 can be calculated 
unambiguously. 

The corresponding diagrams with photons and massless neutrinos in the loop (see Fig. |TJ) 
were calculated by Radkowski ||, Capper, Duff, and Halpern 0, Capper and Duff ||, Duff 
and Liu ||. This contribution to the numerical factor a qu is 



3 Still, erroneous papers on the subject are being published up to now. We mean the articles ||, where the 
claim is made that the classical relativistic corrections to the equations of motion of two bodies separated by 
large distance depend essentially on the inner structure of these bodies. We believe that this claim does not 
withstand criticisms. 




Figure 1: Photon (neutrino) loop 



where N v is the number of massless two-component neutrinos. 

As to the contribution to the effect from the graviton exchange, it was considered by 
Donoghue [10-13], Muzinich and Vokos [14], Hamber and Liu [[0|], Akhundov, Belucci and 



Shiekh |16| |. However, there are no quantitative agreement among the results of these works, 
even the predictions for the sign of the correction differ. 

We believe that the correct result for the quantum correction to the Newton law will be 
sufficiently interesting from the theoretical point of view. This is the aim of our investigation. 
Among the previous works on the subject, the most detailed presentation of the calculation is 



given in [11], [Lq| . Our approach — the direct calulation of Feynman diagrams, the choice of 
the field operator for the gravitational field and of the gauge — is the same as in [10-13, 16]. 
It allows for a relatively detailed comparison of calculations of separate contributions to the 
effect. This comparison has demonstrated that in [10-13, 16] not all diagrams are taken into 
account, and the considered contributions are calculated incorrectly. Below, when discussing 
concrete diagrams, we will come back to the comparison with the previous works, including 



14] , |15 . And meanwhile, let us note an obvious error in [10-13, 16]: therein formula for the 
Fourier-transform of the function 1/r 3 (see (^)) contains n 2 , instead of n, and this error persists 
in the final answer as well. 

Some of the considered diagrams contribute also to the classical relativistic correction. To 
check our calculations we computed in parallel these classical contributions and compared 
them with the corresponding results of ||. As to these classical corrections, we have complete 
agreement with M for each diagram taken separately. 

2 Propagators and vertices 

We use below the units with c = 1, K — 1. 

As a field operator h^ v we choose the deviation of the metrics from the flat one: 

9im> = 5fiv + xh^ u ] 5^ = diag(l, -1, -1, -1); x 2 = 32irk = 32nl 2 . (5) 



We use the gauge where the graviton propagator is 

P B 1 

The tensor P^^p is conveniently presented as [ |Hj l 

where I^ U)a p = \ (^a^/3 + ^0^/3) is a sort of a unit operator with the property 

1 T jxv ,a{)t 'a/3 ^fiu 
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Figure 2: Vertices 



for any symmetric tensor t a p. Let us note the following useful identity: 

Propagator of a scalar particle is the usual one: 

G(p) ' 



p 2 + i0 



Vertex of the interaction of a scalar particle with graviton is (see Fig. |2]a) 



m 



The contact interaction of a scalar particle with two gravitons (see Fig. [|b) is 



V>cX,pa = il* 2 [ Ix\,a5hp,pa{PaPp + PaPp) ~ ^ (^xIpa,a/3 + S pa I>eX,ap)PaP'p 
(P' ~ P? 



+ 
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pa 



(7) 
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(9) 



(10) 



To our accuracy, one can neglect in this expression the last term, with (p' — p) 2 . 

Let us note that in the works [IT, |TB[ the vertex (PH|) is erroneously presented (and indeed 
used in the calculations) with a factor two times smaller, x 2 /2 instead of x 2 . We will come 
back to this factor in Section 4. 

The following useful identities are worth mentioning here: 



Ppu, a p [PaP'p + p' a Pf3 ~ 5 a p(pp' - m 2 )] = p a p' p + p' a pp - 5 Q/9 m 2 

PatP,*cxP'y5,pcr**cX,p(T • 



(11) 

(12) 



As to the 3-graviton vertex (see Fig. 0c), which has the most complicated form, we will 

(13) 



follow and present it as 

T/ - • — ' 



'VoW = -Pa/3,75 [kpK + (k — q)n(k - q) v + q^q v - - S^q 2 ], 
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Figure 3: Simple loops 
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In this vertex as well one can neglect, to our accuracy, the last structure 5 v fJ , 1/ja p jr s- 

Let us note that in papers fTl], [16|] vertex ([H|) is written down with an opposite sign. Our 
sign is confirmed by the following argument: for physical gravitons with momenta k and k — q 
vertex fllTf ) should agree in the limit q — > with the interaction @ of a graviton with a scalar 
particle. 



3 Simple loops 

It is convenient to start with the diagrams where the Feynman integrals contain two denomi- 
nators only. 

The simplest of them, Fig. [3]a, is missed at all in [10-13, 16]. Its calculation causes no 
difficulty, one has only to use the identity (|P2|) and to go over to the nonrelativistic limit in 
both 2-graviton vertices. The result for this contribution to the quantum correction is 

7T r 6 

The calculation of the next diagram, Fig. ^jp, and that obtained from it by interchanging 
scalar particles, is also sufficiently simple and results in 



26 Pmimi 

U qu2 =- (15) 

Sir r A 



The result of [11] for this contribution differs from (|T^) only by a wrong power of tt. The 



corresponding result of |i| is quite different. 

As to the diagrams Fig. 0c, d with the polarization operator of graviton, we do not have 
much to add to works [10-13] on the method of the calculation, and nothing to add at all on 



the result itself (which is used also in [[U|, |16||). However, for the completeness we present here 
briefly this calculation. 





a b 
Figure 4: Triangle diagrams 



The effective Lagrangian corresponding to the sum of these diagrams, with gravitons and 
vector ghosts, as obtained by 't Hooft and Veltman [[L7fl, is 



L= l —\Yi\q 2 \(—R 2 + —R„ v R^\. (16) 

16tt 2 |y 1 V120 20 M J y 1 

To linear approximation the Ricci tensor and scalar curvature of the external gravitational field 
which enter this expression are 

R = xh a p(q 2 5 a( i - q a qp) = xh a/3 r al3 . 
The corresponding contribution to the graviton polarization operator is 

n Q/3 , 7 5 = - ^ In I q 2 \ + ^ r^apr^A . (17) 

We have taken into account here two possibilities of identifying R^ and R with the upper and 
lower external gravitons. The subsequent calculation is straightforward. Let us mention only 
that the summation over fi, v is conveniently performed at the end. Finally, this contribution 
to the quantum correction is 

43 k 2 mi m 2 

Uqu3 = - ™ § • (lo) 

30tt r 6 

Let us mention that diagrams ^c,d were computed in other variables, ip^ = h^ u — l/25^ u h, 



in Refs. |Tq, 18, 19] , and in the source description of gravity due to Schwinger in 



4 Triangle diagrams 

A sort of master formula for triangle diagrams, Figs. f|a,b, reads (we keep here only terms 
singular in \q\): 

d^k 1 1 i 7r 2 m , . ., . . 

+ In I q I • (19) 



(2tt) 4 k 2 (k - q) 2 ((p - k) 2 - m 2 ) 32ir 2 m 2 \ 

It is conveniently obtained by calculating first the imaginary part of its lhs in the t-channel, 
and then restoring its rhs via the dispersion relation. The first term in the rhs of formula 
( [HID generates 1/r 2 in the coordinate representation and contributes to the classical relativistic 



Figure 5: Box diagrams 



correction. It is kept in (|19|) to check the calculations by comparison with the corresponding 
results of j§. 

Our result for the contribution of more simple diagrams of the type Fig. |]a is 

28 k 2 mimo 

U quA = ^ . 20 

7r r 6 



This contribution is also missed in [11, [16|] 



These diagrams contribute to the classical correction as well. An extra proof of our normal- 
ization for the seagull vertex is the agreement with the corresponding classical result of 0. 
Much more tedious is the calculation of diagrams of the type Fig. |]b. It results in 

_ 29 fc 2 mim 2 



The corresponding result of [|16[ differs from fl2~T|) only by sign and wrong power of 7r. The result 
of |TTJ for this contribution is quite different. 

5 Box diagrams 

The expressions for the matrix elements corresponding to the box diagrams Figs. |5]a,b can be 
presented as 

M =i ^ [ d ' k [a-b{D^D,)}\ 

16mim 2 J (2tt) 4 k 2 (k - q) 2 D x D 2 ' V ; 

D 1 = k 2 -2( Pl k), D 2 = k 2 + 2p 2 k, a = 2( PlP2 ) 2 - m 4 , b=( Pl p 2 ). 



r d 4 k [a' + b'(D 1 + D' 2 ) } 2 

16mim 2 J (2vr) 4 k 2 {k - q) 2 D x D' 2 



D' 2 = k 2 -2(p 2 + q,k), b'= (pi,p 2 + q), a' = 2( Pl ,p 2 + q) 2 - m 4 . 

It is convenient to single out in the numerators of these integrals such structures that cancel 
one or both denominators D x , D 2 (D' 2 ). When cancelling a single denominator, one is left with 
effectively triangle diagrams of the type Fig. [|a: 



A 



d 4 k (. 2 f2(j, 2 k) 2( Pl k)\ I 1 1 



Mi " = 'm^J wm^wv \~ — w) +2ab \Di + vj\- (25) 

It can be easily demonstrated that to the accuracy we are interested in, expressions (p3|) and 
( f2~5| ) cancel. 

Expressions with two denominators deleted, 

x 4 f d 4 k 2b 2 
2s ~ 1 16m im2 J (2^ k\k-q) 2 ' (26) 

2u ~ 1 16m im2 J (2^0* k 2 {k-q) 2 ' ( } 

correspond to the diagrams of the type Fig. ||a. These contributions add up into the following 
result for the effect discussed: 

t^-**^. (28) 

Now we are left with the "irreducible" parts of diagrams Figs. ||a,b. These irreducible 
matrix elements are conveniently obtained by calculating first their imaginary parts, in the s 
and u channels respectively, and then restoring the real parts through the dispersion relations. 
The results are (we omit nonsingular in | q 2 | terms) 

. x 4 a 2 r d A k 1 
0s ~ 1 lQ mi m 2 J (2^ (k 2 - \ 2 ){{k - qf - A 2 )(P - 2 Pl k){k 2 + 2p 2 k) ^ ' 
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In ^ ; (30) 



. x 4 a 2 r d 4 k 1 
0u ~ % 16m im2 J J2^Y {k 2 - \ 2 ){{k - q) 2 - \ 2 ){k 2 - 2 Pl k){k 2 - 2{p 2 + q, k) 

\q 2 \ 

lniy. (32) 
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In the above formulae s = (pi+p 2 ) 2 and u = (pi — p 2 — q) 2 . Expressions (|2"5|), ( [5T] ) are conver- 
gent in the ultraviolet sense, but diverge in the infrared limit, depending logarithmically on the 
"graviton mass" A. As usual, such behaviour is directly related to the necessity to cancel the in- 
frared divergence in the Bremsstrahlung diagrams (of course, the gravitational Bremsstrahlung 
in the present case). Though the leading singularity in q is of the type In | q 2 \/\ q 2 \, a term with 
In | q 2 \ arises in the sum of the irreducible boxes as well, and generates the following quantum 
correction to the Newton potential: 



23 k 2 m 1 m 2 

Uqu7 — - t. n — • (odj 

Sir r 6 



It is worth mentioning that, as distinct from the previous contributions where \q 2 \ served as 
an infrared cut-off for ultraviolet divergent integrals, here | q 2 \ is the upper limit for infrared 
divergent integrals. 

For the box diagrams as well, we have checked that our results for thus generated classical 
corrections agree completely with those of [|J. 

The box contributions to the quantum correction are missed at all in [10-13, 16], though 
diagrams Fig. |^a,b are considered in [M from another point of view. 



On the other hand, neither in |L4], nor in [15] we could find any mention of the "infrared" 
contribution of the type ([!3|). In fact, in [|TJ] the problem of classical and quantum corrections 
was treated in different variables, = ~ 1/25^/1. It can be easily demonstrated that the 
expressions for the box diagrams are exactly the same in both variables, ip and h. However the 

]15| , disagree both with the classical ones obtained in j| 
to be the same in both variables, ip and h) and with 
or fl33[), or the sum of (EBT) and (1331). 



box contributions, as calculated in 
(which are demonstrated explicitly in 
our results for the quantum correction, be it 

At last few words more on Ref. [14|. The approach advocated therein looks quite interesting 



and promising. However, the results for the quantum correction presented in [|TJ] do not agree 
with ours (neither do they agree with those of [10-13, 15, 16]). Due to the lack of details in 
we cannot say with certainty what is the origin of the disagreement. Still, an impression 



14 



arises that at least it is overlooked in [14| that the irreducible triangle diagrams generate not 
only classical corrections, but quantum corrections as well, i.e. it seems that in [H] the second 
term is omitted in formula (|19|). 



6 Conclusions 



Summing up all the contributions obtained, (]14l) , (|15l) , (|18l) , (|20D, (pij), fl28|), (|33|), we arrive 
at the following result for the Newton potential, with the discussed quantum correction due to 
the two-graviton exchange included: 



U{r) 



1 + 



121 kh 



107T c 3 r 2 



(34) 



Let us note that the derived overall correction enhances, but not suppresses the common Newton 
attraction. 
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